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An abrupt change in a decay systematics around the A*' — 126 neutron shell closure is discussed. 
It is explained as a sudden hindrance of the clustering of the nucleons that eventually form the a 
particle. This is because the clustering induced by the pairing mode acting upon the four nucleons is 
inhibited if the configuration space does not allow a proper manifestation of the pairing collectivity. 
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I. INTRODUCTION 

It is nearly a century ago that the Geiger-Nuttall law, 
which was to revolutionize physics by its implications, 
was formulated based on a decay systematics [l], Q ■ In- 
deed, its explanation by Gamow Q and also by Gurney 
and Condon [1] required to accept the probabilistic inter- 
pretation of Quantum Mechanics. The extend to which 
this was revolutionary can perhaps best be gauged by 
noticing the multitude of models that have been put for- 
ward by outstanding physicists as an alternative to the 
probabilistic interpretation. This debate rages even at 
present 

The Gamow theory reproduced the Geiger-Nuttall law 
nicely. One can assert that this is an effective theory, 
where concepts like "frequency of escape attempts" have 
to be introduced. Yet Gamow's theory is so successful 
that even today it is applied, with minor changes, in the 
studies of radioactive decays (e.g., Refs. In fact, a 

proper calculation of the decay process needs to address 
first the clustering of the nucleons at a certain distance 
outside the nuclear surface and, in a second step, the 
evaluation of the penetrability through the Coulomb and 
centrifugal barriers should be performed at the distance 
where the cluster was formed. The first step is a chal- 
lenging undertaking because a proper description of the 
cluster in terms of its components requires a microscopic 
many-body framework that is very complicated. This 
is the reason why usually effective approaches are used 
when dealing with clusterization. That is, one evaluates 
the penetrability, which is an easy task specially if semi- 
classical approaches are applied, and free parameters are 
introduced for the clustering process trying to reproduce 
experimental data. 

One may then wonder why effective approaches have 
been so successful. The reason is that the a-particle for- 
mation probability usually varies from nucleus to nucleus 
much less than the penetrability. In the logarithm scale 
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of the Geiger-Nuttall law the differences in the forma- 
tion probabilities are usually small fiuctuations along the 
straight lines predicted by that law Q for different iso- 
topic chains. The importance of a proper treatment of a 
decay was attested by a recent calculation which shows 
that the different lines can be merged in a single line. One 
thus obtained a generalization of the Geiger-Nuttall law 
which holds for all isotopic chains and all cluster radioac- 
tivities [l^ [ll|. In this universal decay law (UDL) the 
penetrability is still a dominant quantity. By using three 
free parameters only, one finds that all known ground- 
state to ground-state radioactive decays are explained 
rather well. This good agreement is a consequence of 
the smooth transition in the nuclear structure that is of- 
ten found when going from a nucleus to its neighboring 
nuclei. This is also the reason why, e.g., the BCS approx- 
imation works so well in many nuclear regions. 

In this paper we will show that, when a sudden tran- 
sition occurs in a given chain of nuclei, departures from 
the UDL can be seen. Our aim is to understand why this 
difference appears. We will also try to discern whether 
one can. in general, obtain information about the struc- 
ture of the nuclei involved in the decay. This would be 
an important task because many regions of the nuclidic 
chart now under scrutiny, especially superheavy nuclei, 
are radioactive and often a decay is the only tool that 
one has to explore their structure. 

In Section |ll] the formation amplitude is defined. In 
Section IIIII a formation amplitudes extracted from ex- 
perimental data are presented and abrupt changes are 
noted. In Section ITVl the evaluation of the formation am- 
plitudes and half-lives of Po isotopes, which do not follow 
the UDL, is performed. A summary and conclusions are 
in Section IVl 



II. THE FORMATION AMPLITUDE 

After the seminal Gamow's paper, the first attempt to 
formulate a proper treatment of a decay was based on the 
compound system theory developed by Teichmann and 
Wigner [l^ . Here the very complicated process occur- 
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ring as the compound system decays, is divided into an 
"internal region" , where the compound state is restricted, 
and the complementary "external region" . This division 
is such that in the external region only the Coulomb and 
centrifugal forces are important. Thus the decaying sys- 
tem behaves like a two-particle system. This formulation 
was applied by Thomas to a-decay [Tsj to obtain the 
classical expression for the decay width F/ as 



Vi{R) = 2Vi{R) — \Fi{R)\\ 



(1) 



regular Coulomb functions, respectively. From Eq. ([T]) 
it is straightforward to see that the width Fj [R) cannot 
depend upon R, since outside the range of the nuclear 
interaction (i.e., just outside the nuclear surface) the in- 
ternal and external wave functions are the same [21|, i.e., 
J-i{R) cx GL{R)+iFi{R). This is of course valid provided 
that the formation amplitude was evaluated properly. In 
fact a way of probing the calculation is just by investi- 
gating whether the width is dependent upon i?,, and in 
such a case by how much (2^ . 

The a-decay half-life can be written as 



where I is the angular momentum carried by the outgo- 
ing a-particle, V is the penetration probability and ^ is 
the reduced mass corresponding to the final system con- 
sisting of an a particle and a daughter nucleus. R is the 
radius dividing the internal and external regions. At this 
point the wave function of the a-particle already formed 
in the internal region is matched with the corresponding 
outgoing two-body wave function in the external region. 
The amplitude of the wave function in the internal region 
is the formation amplitude, i.e., 

Fl{R)^ I dRdeddC„[*(Crf)0(C„)y,(R)]r/„M„*m(&,Ca,R) 



(2) 

where d, a and m label the daughter, a particle and 
mother nuclei, respectively. 4* are the intrinsic wave 
functions and ^ the corresponding intrinsic coordinates. 
4'(£,a) is a Gaussian function of the relative coordinates 
of the two neutrons and two protons that constitute the 
a-particle, coupled to zero angular momentum [isj . 
The rest of the notation is standard. 

One sees from Eq. ([2|) that [R) would indeed be the 
wave function of the outgoing a particle ipa{R) if the 
mother nucleus would behave at the point R as 



(3) 



Since this is usually a smal component of the mother 
nucleus wave function, the corresponding formation am- 
plitude ^ is small, of the order of lO'^ [Hi. The mam 
problem in the evaluation of this quantity is the descrip- 
tion of the clusterization of the four nucleons that even- 
tually become the a-particle. In pursuing this task one 
has found that the mode that determines clusterization 
is the pairing vibration [l^ In fact, the study of 

a-clusterization gave rise to the realization that there 
should be a giant pairing vibration lying high in the nu- 
clear spectra [3 113 • It is also interesting to notice that 
the a-clusterization in a-decaying nuclei has triggered 
the appearance of effective models where the wave func- 
tion of nuclei such as ^^^Po is assumed to have the form 
([3]). The spectra thus obtaine d ag ree well with the cor- 
responding experimental data [20| . 

Going back to Eq. ([1]), the wave function correspond- 
ing to the external region, i.e., to the outgoing chan- 
nel, gives rise to the penetration probability ViiR) = 
kR/{Gf + F^), where Gi and Fi arc the irregular and 
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where v is the outgoing velocity of the emitted particle. 
The distance R will be taken as the touching point, i.e., 
R = Ro{A]^^ + 4^/3), with i?o=1.2 fm. The other quan- 
tities arc standard, i.e., H^{x,p) is the Coulomb-Hankel 
function with arguments x ~ 'iZde'^/hv and p = p.vR/h. 

In microscopic theories the formation amplitude is 
evaluated starting from the single-particle degrees of free- 
dom of the neutrons and protons that eventually become 
, the cluster. This requires advanced computing facilities 
as well as suitable theoretical schemes to describe the 
clustering process. It is therefore not surprising that the 
first calculations of absolute decay widths were performed 
after the appearance of the shell model. These calcula- 
tions had limited success due to the small shell model 
spaces that could be included at that time [ij]. Yet, in 
retrospect it is surprising to note the deep insight the pi- 
oneers in these shell model calculations had on the role 
of configuration mixing to induce clustering [23j . That 
this was indeed the case was shown much later [l^ 
in the case of the decay of the nucleus ^^^Po with two 
protons and two neutrons outside the doubly magic core, 
^"^Pb, which has been considered as a textbook exam- 
ple in illustrating the clustering and decayof the alpha 
particle in heavy nuclei (see, e.g., Ref. [2J| and refer- 
ences therein). In fact this case is very important for the 
present paper, since the most significant departure of the 
UDL from experimental data that we will investigate is in 
the ground-state to ground-state decays of Po isotopes. 



III. THE EXPERIMENTAL DATA 

Using the experimental decay half-lives [25| one can 
extract the formations amplitudes by applying Eq. 
One thus obtains 
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This is shown in Fig. [T] for different even-even isotopes 

as a function of the quantity p' ~ \J 2AZd{A]/'^ + 4^/^), 
where A = AAd/Am- This is one of the two variables 
that defines the UDL [l^. 
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FIG. 1. (Color online) log^j \RF{R)f as a function of p' . The 
solid line denotes the smooth behavior of the UDL. The values 
between the two dashed lines differ from the corresponding 
UDL values by a factor of at most three. 



One notices in Fig. [T] that at p' « 70 a division occurs 
between decays corresponding to iV < 126 and N > 126. 
Perhaps even more important is that for most cases the 
UDL predicts the experimental values within a factor of 
three, except for N = 126, where the difference becomes 
about one order of magnitude. This is so distinct that 
one may even suspect that the difference in the values of 
J^a when going from one nucleus to its neighbors in the 
vicinity of N = 126 overruns the corresponding differ- 
ences in the penetrability. If one understands the reason 
of this large variation, a decay may provide a powerful 
tool to study the structure of decaying nuclei. This point 
will be analyzed in the next Section. 

The case that shows the most significant hindrance 
corresponds to the a decay of the nucleus ^^°Po, 
with logio I i?J"a(i?)P < -3 fm"\ The symbols with 
logiQ \RTa{R)\^ ^ —2.7 fm^^ correspond to the a de- 
cays of nuclei ^°^Po [N = 124), ^^^Rn {N = 126) and 

194pb. 

Finally, it is worthwhile to point out that this sudden 
change in ct-deca y sy stematics at TV = 126 has also been 
noticed in Refs. [26l - l29i . Moreover, in the semiclassical 
approaches of Refs. [30M36l| the decay half-lives of nuclei 
with N = 126 are significantly underestimated. In the 
alpha-decay formula of Refs. [3^ [35j . an empirical cor- 
rection term has been introduced to taken into account 
the large underestimation around shell closures. 



IV. EVALUATION OF THE HALF-LIVES IN PO 
ISOTOPES. 

In this Section we will analyze, within a microscopic 
formalism, the half-lives of the isotopes that show the 
kink at = 126 discussed above. Wc will take the de- 



cay of ^^°Po as a typical example and compare it with 
that of ^^^Po. To compare with experimental data, we 
extract the magnitude of the formation amplitudes from 
measured half- lives by using Eq. ([S]). One thus obtains 
the value J'aiR) = 3.305 x IQ-^ fm-3/2 in ^lOpo and 

J"a(i?) = 1.082 X 10~2 fm-3/2 212pQ^ ^^^^.^ ^g^j 

i? = 9.0 fm. These correspond to a variation in the for- 
mation amplitudes by a factor of 3.28, that is a factor of 
10.73 in the formation probabilities. 

Within the shell model a four-particle state 04 in ^^^Po 
can be written as 

|2i2po(a4)) = J2 ^(«2/32;a4)P'°Pb(a2)®''°Po(/32)> 

(6) 

where a2 {132) labels two-neutron (two-proton) states. 
For the ground state of ^^^Po, it was found (2^ 
that X(2iopb(gs) (E) 2iopo(gg)) ^ 0.9, while 
X(2iopb(2+) (g) 2iopo(2+)) ^ _o.3. 

Each of the terms in Eq. ([6]) corresponds to neutron- 
neutron (nn) or proton-proton (pp) states, i.e., states 
determined by the nn or pp interaction. The neutron- 
proton (up) interaction mixes those states. In other 
words, the amplitudes X are influenced by the np in- 
teraction. If this interaction is neglected, then only one 
of the configurations in Eq. ^ would appear. This is 
done, for instance, in cases where the correlated four- 
particle state is assumed to be provided by collective vi- 
brational states. Rather typical examples of such states 
are p^°Pb(gs)) and p^°Po(gs)). It is therefore not sur- 
prising that calculations have been performed by assum- 
ing that p^2Po(gs)) is a double pairing vibration (Tsl-ISTt. 



'Po(gs)) 



'Pb(gs) . 



'Po(gs)). 



(7) 



The corresponding formation amplitude acquires the 
form, 

J-„(i?;2i2po(gs)) = J dKd^MU 

X vl/(ri, r2;2i" Pb(gs))*(r3, r4;2i° Po(gs)), (8) 

where ri, r2 (rs, r4) are the neutron (proton) coordinates 
and R is the center of mass of the a particle. 

With this expression for the formation amplitude the 
experimental half-life is reproduced rather well if a large 
number of high-lying configurations is included. These 
configurations are needed to describe the clusterization 
between the two neutrons and the two protons in the a 
particle. Yet the corresponding a-dccay half-life is still 
too small by more than one order of magnitude. This is 
because the neutron-proton interaction is not included in 
Eq. ([7]). When this is done, and again a large configu- 
rations space is used, the neutrons and protons become 
also clustered, enhancing the value of the half-life. It is 
also important to underline that the inclusion of the large 
configuration space provides a half-life which is indepen- 
dent upon the matching point R [2^ . 
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We reproduced these calculations by using a surface 
delta interaction and nine major shells of a harmonic 
oscillator (HO) representation. The decay of the nucleus 
^^°Po(gs) leads to the daughter nucleus ^°^Pb(gs), which 
is a two-hole state. Here we used the five HO major shells 
corresponding to the single-hole states that describe the 
wave function of ^°^Pb(gs) as 



1 206 



Pb(gs)) = ^ X(/ii/l2;206pb(gs)) 

hi </i2 

(^^1^/^2)0+ ,208t 



V2 



^Pb(gs)), 



(9) 



where h labels single-hole states and the hole creation 
operator is standard, i.e., 6^ = (— l)-'^'"cj_„i. The for- 
mation amplitude becomes, 

J-„(i?;2iOpo(gs)) = JdRd^MU 

x**(ri, r2;'°^ Pb(gs))(r3, r4;2i" Po(gs)). (10) 



By comparing Eqs. ([8]) and ([TOl) one sees that the 
only difference between the two expressions is the two- 
neutron wave function, which corresponds to the two- 
particle state ^^''Pb(gs) in Eq. ^ and to the two-hole 
state ^°^Pb(gs) in Eq. pUj) . Therefore the kink observed 
experimentally should be related to the difference in clus- 
terization induced by the pairing force in these two cases. 
To analyze the clustering features we will consider only 
the spin-singlet component, i.e., (xiX2)oi of the two-body 
wave function, since that is the only part entering the in- 
trinsic a-particle wave function. This component has the 
form. 



X(pp{ri)(pq{r2)Pipicoa 6*12), 



(11) 



where ip is the single-particle wave function and Pi is the 
Legendrc polynomial of order / satisfying Pi(cosO) = 1 
(notice that for the ground states studied here it is 
Ip = Iq)- As mentioned above, the pairing vibrations 
show strong clustering features as the number of single- 
particle states is increased But another manifesta- 
tion of the pairing collectivity is an enhancement of the 
wave function on the nuclear surface. The reason of this 
enhancement is that all configurations contribute with 
the same phase in the building up of the two-particle 
wave function on the nuclear surface. The same mecha- 
nism increases the a formation amplitude and, therefore, 
the relative values of the wave functions of ^^°Pb(gs), 
^^°Po(gs) and ^°^Pb(gs) on the nuclear surface give a 
measure of the importance of the corresponding forma- 
tion amplitudes. 

To study the behavior of the two-particle wave func- 
tions we will apply Eq. (im with ri = r2 and 6*12 = 0. 
This is reasonable since due to clustering the wave func- 
tion is strongly peaked at 612 = 0. Calling R = ri, we 



have plotted in Fig.[2]^'2(i?, R, 0) as a function of R. One 
sees that the wave functions are indeed strongly enhanced 
at the nuclear surface, as expected. But the important 
feature for us is that the enhancement is strongest in 
^^°Pb(gs) and weakest in ^*'^Pb(gs). This is because there 
is a relatively small number of configurations in the hole- 
hole case. In addition, the radial wave functions corre- 
sponding to the high-lying particle states extend farther 
out in space with respect to the hole configurations. 
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FIG. 2. (Color online) The two-body wave function 
'I'2(R,R,0) corresponding to the pairing vibrations in the two- 
neutron particle ^^*^Pb(gs), two-proton particle ^^'^Po(gs) and 
two-neutron hole ^''^Pb(gs) cases. 

With these two-body wave functions we proceeded to 
evaluate the a formation amplitudes in ^^^Po(gs) and 
^^°Po(gs). The results are shown in Fig. |31 One finds 
from this figure that with i? = 9 fm the observed ratio 
between the formation amplitudes in ^^^Po and ^^°Po can 
be reproduced nicely. 



A. The neutron-proton interaction 

We have assumed (Eq. ([7])) that ^^^Po(gs) is virtually 
a correlated two-neutron two-proton state. The same is 
valid for ^^°Po(gs), although here the state is a correlated 
two-particle (proton) two-hole (neutron) state. This is a 
manifestation of the pairing vibrational character of two- 
particle states in the Pb region. That is, the correlated 
two-particle and two-hole states in the Pb region can be 
considered as boson degrees of freedom. This was one of 
the main assumptions in the Nuclear Field Theory [38| as 
well as in the original Interacting Boson Model |3y| . This 
assumption implies that the neutron-proton interaction 
does not play a very important role in the spectroscopy 
of the states. However, as we have seen, this interac- 
tion induces the clusterization of neutrons and protons. 
As pointed out in Ref. [l^, in the Pb region low lying 
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V. SUMMARY AND CONCLUSIONS 
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FIG. 3. (Color online) The a formation amplitudes RTa{R) 
corresponding to the nuclei ^^^Po(gs) and ^^°Po(gs). 



neutron and proton single-particle states are very differ- 
ent from each other, or are particle-hole states. There- 
fore the neutron-proton interaction affects only slightly 
the ground states and the clusterization occurs through 
high-lying configurations. This point is supported by our 
shell-model calculations with the surface delta as well as 
realistic interactions. 

Only when neutrons and protons move in the same 
orbits it is expected that the neutron-proton interaction 
would affect significantly the spectroscopic properties as 
well as the clusterization. We confirmed this by studying 
a model case in which the core consists of an equal num- 
ber of neutrons and protons, namely the a decay of the 
fictional nucleus ^^^Po(gs), with two neutrons and two 
protons outside the core ^^'*Pb(gs). We used, for neu- 
trons as well as protons, the single-particle states cor- 
responding to protons in the study performed above for 
^^°Po. We also used the same interaction. As expected, 
we again found that neutrons and protons are strongly 
clustered as a result of the corresponding pairing inter- 
action. But also the proton-neutron clustering is signifi- 
cantly enhanced by the proton-neutron interaction. This 
indicates that in realistic N = Z nuclear regions, for in- 
stance around ^"^Sn, there should be a large probability 
to form an alpha particle (see, e.g., Ref. [40|). One can 
thus conclude that alpha decay probes may be a power- 
ful tool to get information about the structure of heavy 
N Ki Z nuclei which, otherwise, would be difficult to 
reach. 



In this paper we have applied the recently proposed 
universal decay law (UDL) [l^l to perform a systematic 
calculation of a decay half-lives over all experimentally 
known cases. We found that although the UDL repro- 
duces nicely most available experimental data, as ex- 
pected, there is a case where it fails by a large factor. 
This corresponds to the a decays of nuclei with neutron 
numbers equal to or just below N = 126. The reason for 
this large discrepancy is that in N < 126 nuclei the a 
formation amplitudes are much smaller than the average 
quantity predicted by the UDL (Fig. [T]). This is an in- 
dication that the a decay transitions in these nuclei are 
hindered with respect to those in the open-shell region. 

The case that shows the most significant hindrance cor- 
responds to the a decay of the nucleus ^^"Po for which 
standard shell-model calculation is feasible. Starting 
from the formal definition of Eq. 1^ , we calculated the a 
formation amplitude of ^^"Po and compared it with that 
of ^^^Po. In these two cases the formation amplitudes 
can be described by the simple expressions ^ and (|10p . 
We found that the formation amplitude in ^^°Po is hin- 
dered with respect to the one in ^^^Po due to the hole 
character of the neutron states in the first case. This is 
a manifestation of the mechanism that induces cluster- 
ization, which is favored by the presence of high-lying 
configurations. Such configurations are more accessible 
in the neutron-particle case of ^^^Po than in the neutron- 
hole case of ^^°Po. This is a general feature in nuclei 
where neutrons and protons occupy different low-lying 
major shells. If instead both types of particles occupy 
the same shells, the neutron-proton interaction is very ef- 
fective to induce clustering and the formation amplitude 
increases strongly. This was the case in a calculation 
that we performed considering the fictitious N = Z = 
i68pQ igotope as the mother nucleus, indicating that even 
in physically meaningful N = Z nuclear region a decay 
can be enhanced by large factors. 

This allows one to assert that a decay is a powerful tool 
to investigate the shell structure of very unstable nuclei 
(including superheavy ones), where often only a-decay 
quantities can be measured. 
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